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• Abstract. We compute the phase transition curves for mesoscopic supercon- 
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ductors. Special emphasis is given to the limiting shape of the curve when the 
magnetic flux is large. We derive an asymptotic formula for the ground state of 
the Schrodinger equation in the presence of large applied flux. The expansion is 
shown to be sensitive to the smoothness of the domain. The theoretical results 
are compared to recent experiments. 
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In their classical experiment, Little and Parks 0] considered the transition between normal 
^ ! and superconducting phases in cylindrical shells. They discovered that the phase boundary in the 

in ■ 

{Sj ' T — H plane, where T is the temperature and H is the applied magnetic field, is oscillatory. When 

m : 

O ■ the thickness of the shell is very small relative to the coherence length, the period is approximately 
0O ! the fundamental flux quantum O = h/e*c, where e* = 2e is the charge of the Cooper pair. Taking 



into account the effect of the finite thickness of the ring, it can be shown theoretically |TD[ and 
verified experimentally (e.g. [12lJ13|) that the oscillations in the T — H line are superimposed on 
a parabolic background. The situation becomes more complex if the circular symmetry of the shell 



i 

q \ is broken 0, but the T — H line is still oscillatory. 



The classical literature on the subject assumes that the superconducting sample is multiply 



connected, typically with a ring-like geometry. Nevertheless, an oscillatory T — H line is also 
obtained for simply connected domains. This has been shown experimentally by Buisson et al. || 
who considered a small disc. Later Moshchalkov et al. [|1J and Bruyndoncx et al. measured 
the transition temperature T as a function of the applied flux for a mesoscopic square. Both 
groups reported an oscillatory phase boundary superimposed on a linear background. The heuristic 
reasoning for the oscillations is that as the flux increases, the wave function concentrates near 
the boundary, and the sample appears effectively as a thin ring. This gives rise to a topological 
quantization constraint on the phase of the wave function and hence the Little Parks oscillations. 



The analogy with the experiments in thin shells is not perfect, though. In the case of a simply 
connected sample under strong fields, the radius of the 'effective' ring is not fixed. Rather it shrinks 
as the applied field increases. 

We shall show here that the base T — H curve is parabolic for small fluxes, linear for large fluxes, 
and that its precise shape depends on the smoothness of the domain. We compute the T — H curve 
for arbitrary mesoscopic two dimensional samples Q. A constant magnetic field of magnitude H is 
applied in the direction orthogonal to the plane. We pay particular attention to the asymptotic limit 
in which the flux is large. Near the phase transition, the induced magnetic field is negligible and 
the wave function is well-approximated by the solution to the linearized Ginzburg Landau equation 

— {ihV + —A)H = am, (1) 

2m c 

(where V is the two-dimensional gradient operator) with nonreflecting boundary conditions at dfl 
the boundary of Q. Here A is the magnetic vector potential satisfying VxA = Hz, m is the 
effective electron mass, a = -^-j? T ~ Tc , where Co is the coherence length at T = and T c is the 
critical temperature in the absence of magnetic field. Since we are interested in finding the relation 
between T and H, the problem is mathematically equivalent to computing the ground state of the 
Schrodinger equation for a particle under a magnetic field. The strong field asymptotics can also 
be shown to be equivalent to evaluating the critical magnetic field H c3 for materials with large 
Ginzburg Landau parameter 0, @. 

We denote by the applied flux measured in units of 4>o- When cf> ^> 1 we expect \^\ to be very 
small except at a boundary layer whose nondimensional thickness is O(0 - ^). The detailed structure 
of and the asymptotic expansion of a depend upon the geometry and smoothness of Q. 

It is convenient to scale the domain Q canonically through a similarity transformation that 
preserves its shape while changing its area to 2tt. The new domain will be denoted by Q so that we 
must solve the non-dimensionalized problem 

(zV + (f)A N f^ = //* in tt, 

where A^ satisfies V x A^ = lz. When the domain is smooth (i.e. has no corners), the leading 
order term in the expansion for T is linear in the flux: Tc ~ T ~ J^o0- Here R = (|fl|/27r)5, and /jq 



is determined by the eigenvalue problem 

Vv, + (/io - (V ~ Vo /2 ) 2 )*P = 0, ^(0) = 0, Jim = 0, V(0) = 1- (2) 
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This equation was introduced by Saint- James and deGennes [0 in their work on nucleation of 
superconductivity near a half plane. 

Higher order terms in the expansion are computed through a careful boundary layer analysis. 
In general the wave function concentrates near the point of the boundary where the curvature k is 
maximal. It decays exponentially fast away from this point both tangentially along the boundary of 
Vt (on a lengthscale of 0(4>~ lj/8 R)) and radially into the interior of f2 (on a lengthscale of 0{4>~ 1 ^ 2 R)) 
[|J. The expansion for fi takes the form 

T -^Y L = §M> + /^ 1/2 + /^ 1/4 + •••)• (3) 

Here \i\ = —j^ and \ii = ( ^ Kss ^ ^ ; where I = / °° ip 2 and K and k ss are evaluated at the point 
of maximum curvature (with k ss assumed to be strictly negative there). Solving (g) numerically we 
obtain /i = 0.59 and J = 1.312. 

In the special case where Q is a disc, the wave function has radial symmetry. In fact, it is possible 
to obtain an explicit solution to the linearized problem in terms of Kummer functions || . However 
it is still beneficial to write down the asymptotic expansion because it captures simple geometric 
characteristics of the domain. We thus obtain 
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T c -T 



T ^(0.590 + /i 1 0^ + /(0) + ...), (4) 

where f(4>) is an 0(1) oscillatory function of that approaches a periodic function as — > oo. 
We note that a rigorous mathematical treatment of this problem within the context of nonlinear 
bifurcation theory has been carried out in 

When Q is a square, the curvature function is singular and the previous expansion is no longer 
valid. Our computations on the linearized eigenvalue problem indicate that the wave function now 
concentrates near the corners. A simple analysis shows that the leading order term in the expansion 
for T is 

^~|^far0»l. (5) 



The coefficient \i c is the ground state for the Schrodinger operator (zV + An) 2 in the quarter- 
plane. The eigenvalue can be estimated by an upper bound derived from the Rayleigh quotient || . 
Surprisingly the upper bound computed in || is much lower than our numerical computation (see 
below) . 

Expansions of the phase transition curve T(H) for low values of the flux follow from regular 
perturbation theory. Let p be the solution of Poisson equation V 2 p = 1 in Q with homogeneous 
Dirichlet boundary conditions, and set Co = InP- To leading order in we find 



Tc-T S 



For example, we obtain Co = 0.25 for a disc, and Co — 0.22 for a square. 

In Figure 1 we present the numerical computations of the values of ^ T srJ as a function of 6 for 
a square (solid line) and for a disc (arrows) with the same area where <fi lies in the range [0, 16]. The 
results are in perfect agreement with our theory for low 0. It is observed that superconductivity 
nucleates for a square at a higher temperature than for a disc. Our asymptotic results indicate 
that the nucleation in a square will be at higher temperature than for any smooth domain in the 
large flux regime. ^From our numerical calculations we estimate fi c = 0.55. As was mentioned 
above this is about twice the estimate of ||. We cannot explain this large discrepancy. We point 



out, however that Kato et al. |Tl| have solved numercially the time dependent Ginzburg Landau 
equation in a square and studied the existence of superconductivity above H c2 . They expected to 



get the behavior predicted in ||14|| . Instead they found what they called an 'anomalous phenomena', 
which, formulated in our notation, is equivalent 0.5 < /i c < 0.555. This value is in agreement with 
our results. 

We compare our expansions with experimental results. We will make a least-squares fit of the 
experimental data to an expansion of the form 



T c -T g 
T c ~ R 2 ' 
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» 1, (7) 



and 
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T c -T £5 



;/; fM , (8) 
To determine the coefficients fi® and ( e from a given experiment we need to know £o an d the actual 



size of the domain. Therefore it is interesting to note the ratio S = fi^/Ce is independent of both £o 
and the domain scale. It is controlled solely by the shape of the domain. 

We first consider the experiment of Buisson et al. || who studied the T — H line for an aluminum 
disc of radius 7.2 fim. For a disc we obtain the theoretical value of the ratio defined above to be 
5o = .59/Co = 2.36. The experimental value is found to be 5 e = 2.26. 

It has been suspected ||,[[| that some of the anomalies in the experiment were caused by the 
contact leads that modified the spherical geometry. Motivated by the theoretical expansion @ we 
shall test this hypothesis by assuming 

+ W), (9) 

In particular, we wish to estimate the coefficient of the 1 / 2 term in the expansion ([5]) since this term 
depends linearly on the maximal curvature. We used the value £o = 0.2 fim obtained |J from the 
low cf) expansion. We found that the radius of curvature at the nucleation point is 6.5 fim ± 0.8 fim 
which is in reasonable agreement with the actual radius of 7.2 fim. Thus we conclude that at least 
at the high regime the results were not affected by the contact leads. 

Moshchalkov et al. |n| performed a similar experiment with a 1 x 1 fim aluminum square. They 
used fluxes of up to 4 flux quanta. The experimental data was fitted to the numerical solution of 
the associated eigenvalue problem for a disc with area 1. For this purpose they chose £o = 0.1 fim. 
A similar comparison was later performed in 0]. A recent experiment from the same laboratory 
H provides more extensive data. In this experiment a 2 x 2 fim square was used. The authors 
measured the phase transition curve for fluxes of up to 16 flux quanta and again compared the data 
to the expected theoretical values for a disc of comparable area. 

We shall compare the data of || for the square to the asymptotic theoretical formulas ([5]) and 
©. Again this comparison is complicated by the need to estimate the value of £o- As a first 
attempt, we estimate the value of £o by comparing the experimental data for low values of to the 
theoretical expansion (|6|) and obtain £o = 0.095 fim. However, using this value, the experimental 
results H predict fi Q e > 0.81 in the large <j) expansion (|7|). This is way above even the theoretical 
value for the disc. In fact, there is no £o that allows one to accurately match the theoretical curve 
to the experimental data. To demonstrate this point, we note that our theory gives the estimate 
5 = 2.5, while the experimental estimate is 5 = 3.7. (Our estimate of 5 is based upon our numerical 



evidence that \l c from expansion (|5|) is approximately 0.55) We point out that other laboratories 
have recently reported [|I6 |, H], || values for £o for aluminum near 0.2 fim. 

One plausible explanation for the deviation lies in the importance of the nonlocal effects due 



to the electric contacts |TJ|]. In the low regime the coherence length is large compared to the 
size of the square, so that wave function, and thus the transition temperature, are affected by the 



neighborhood of the square. Indeed Figure 3 of JT5] indicates that the electric contacts do affect 
the phase transition curve for low 0, and their effect diminishes as increases. Therefore, we 
shall estimate £o using the high data. We thus obtain £o — 0.1195 fim. The circles in Figure 1 
corresponds to the experimental data of H] with this choice of £o- The nonlocal effect was apparently 
less pronounced in || since their domain was 10 times larger. 

An examination of the three curves in Figure 1 reveals that while the experimental curve lies well 
below the theoretical curves, the location of the peaks in the experimental curve perfectly match the 
location of the peaks of the theoretical curve for the square. On the other hand, the experimental 
peaks deviate considerably from the theoretical peaks for the disc. Bruyndoncx et al. introduced 
the notion of 'effective square area' (which is smaller than the true geometrical area) in order to 
obtain a fit between the experimental data and the theory for the disc. We have shown here that 
the actual experimental peaks are in perfect agreement with the theory for the square. We also 
point out that while the experimental slope for low (near 0.5) is about 0.13 (well below the 
theoretical 0.22 value), the experimental curve steepens just before the first peak (0 near 1.7), and 
the local slope is about 0.19. It still remains to explain the deviation between the experimental and 
theoretical curves for intermediate values of 0, both in actual values and in the amplitude of the 
oscillations. A possible explanation is in the way in which the transition temperature is defined. 
Bruyndoncx et al. || set it to be the temperature for which the resistivity drops to half the normal 
value. As the applied flux increases, the temperature drops, and the uncertainty interval shortens. 

To summarize, we have developed general formula (pj), (|j) and (|5]) for the asymptotic transition 
line in mesoscopic superconducting samples in the limit of large flux. We further derived a formula 
(0) for the asymptotics of that curve in the limit of low flux. In the large flux limit and for smooth 
domains, our formula captures the universal linear base curve that relates to the fact that the wave 
function concentrates near the boundary which appears locally as a half plane. It also contains 



terms proportional to fractional powers of the applied flux that depend on universal parameters 
and on the local geometry of the sample's surface. In the special case of a disc our formula resolves 
an oscillatory dependence of T upon <fi. We have also shown that a square sample exhibits a different 
behavior from smooth domains already at the linear leading order term. Finally we used the new 
theory to analyze recent experiments. 
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